Muon diffusion and electronic magnetism in Y 2 Ti 2 07 
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We report a /iSR study in a Y2Ti207 single crystal. We observe slow local field fluctuations at 
low temperature which become faster as the temperature is increased. Our analysis suggests that 
muon diffusion is present in this system and becomes small below 40 K and therefore incoherent. A 
surprisingly strong electronic magnetic signal is observed with features typical for muons thermally 
diffusing towards magnetic traps below « 100 K and released from them above this temperature. 
We attribute the traps to Ti 3+ defects in the diluted limit. Our observations are highly relevant to 
the persistent spin dynamics debate on i?2Ti207 pyrochlores and their crystal quality. 

PACS numbers: 76.75.+i, 75.10.Jm, 75.40.Gb 



The rare-earth titanates and stanates series of com- 
pounds, R2M2O7 (R is a rare-earth ion and M is Ti 
or Sn), which crystallize in the pyrochlore crystal struc- 
ture (space group Fd3m), are prone to strong geomet- 
rical frustration. pQ Their study has revealed a wealth 
of exotic magnetic properties. These include (i) the 
spin-ice ground state of Ho2Ti 2 07 and Dy 2 Ti 2 07,[2]- 
4 (ii) the ground state reached by Yb2Ti 2 07 after a 
sharp transition in the spin dynamics finger-printed by 
a pronounced peak in the specific heat, [5] (iii) the un- 
conventional dynamical ground state of Tb2Sn 2 07 for 
which magnetic Bragg reflections are observed by neu- 
tron diffraction, 6J while no spontaneous magnetic field 
is measured by the zero field (ZF) muon spin relax- 
ation (/iSR) technique, [7] (iv) the persistent spin dynam- 
ics (PSD) detected in the ordered states of Gd 2 Sn 2 C>7, 
Gd 2 Ti 2 7 and Er 2 Ti 2 7 ,[HHIl] and the non-ordered 
PSD state of Tb9Ti9O7.pi [II] 

The first report of PSD was for SrCrgGa^ig using 
/xSR.[15] It was found that this kagome compound does 
not exhibit magnetic Bragg reflections down to 50 mK. 
Although mostly observed by /iSR, PSD has been also 
proposed by other techniques, e.g., in Gd 2 Sn 2 C>7 using 
155 Gd Mossbauer spectroscopy. [5] Both fiSR and the neu- 
tron spin echo have observed PSD in Tb 2 Sn 2 07[Il[T6] as 
well as Dy 2 Ti 2 07.[4j [17j \TE\ However, the existence and 
nature of PSD is still under debate^ UK- Recently, it 
was suggested that PSD does not exist for Dy 2 Ti 2 07 and 
that the observed relaxation results from coherent muon 
diffusion. |20j In that work it was proposed to perform 
/iSR measurements in Y 2 Ti 2 07 to investigate the possi- 
bility of coherent muon diffusion in a compound where 
frustrated magnetism does not play any role as no 4f- 
magnetism should be present. In this letter we report 
ZF and Longitudinal Field (LF) /iSR measurements on 
Y 2 Ti 2 07. We show that there is no detectable neutral 
muonium formation, and that the positive-muon hopping 
rate at low temperature is much smaller than what was 



proposed for Dy 2 Ti 2 07. 20 Also, we have observed a rel- 
atively strong muon spin relaxation which we attribute 
to a small density of Ti 3+ magnetic defects. Note that 
this system has already been studied by Dunsiger using 
/iSR.[2T] Nevertheless, that study was less extensive than 
the one reported here and did not consider the possibility 
of coherent diffusion at low temperature. When a com- 
parison is possible, our data and that from Dunsiger are 
similar. 

Poly crystalline Y 2 Ti 2 07 was prepared by a solid state 
reaction. Starting materials of Y 2 03 and Ti0 2 with 
99.99% purity were mixed and ground. They were then 
heat treated at 900-1150°C in air for more than lOOh 
with several intermediate grindings. The resulting pow- 
der was hydrostatically pressed in the form of rods (8 
mm in diameter and 60 mm in length). The rods were 
subsequently sintered at 1150°C during 15h. The crys- 
tal growth was done using an optical floating zone fur- 
nace with four 1000W halogen lamps as a heat source. 
The growing conditions were: growth rate of 10 mm/h, 
feeding and seeding rods were rotated at about 20 rpm 
in opposite directions (to have homogeneity of the liq- 
uid), growth done in 4 bar pressure of an argon and oxy- 
gen mixture (50:50). The crystal was post-annealed at 
1150°C in argon for 15h in order to remove possible over 
stoichiometric oxygen. Phase purity of the grown crystal 
was checked with conventional powder x-ray diffractome- 
tcr and the obtained lattice parameter a=10.099A is in 
good agreement with the literature. [22] The crystal was 
then aligned using an X-ray Laue camera, and magneti- 
zation measurements were performed down to 1.9 K (see 
Fig.0. 

The ^(SR measurements were carried out on the Gen- 
eral Purpose Spectrometer (GPS) at the Swiss Muon 
Source facility of the Paul Scherrer Institut (Switzer- 
land). Most of the measurements were done in the LF 
geometry for which the initial muon spin and the ex- 
ternal magnetic field B cxt are parallel. [2"5I |2~4"] We define 
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FIG. 1: Inverse susceptibility in SI units as a function of 
temperature in zero field cooled (ZFC) and field cooled (FC) 
protocols. It was verified that the applied field was small 
enough to be in the linear regime. We fit the data to 
Xdc = C CW /(T - Tew) and obtain C C w= 2.45(5) xI(T 4 K 
and Tcw=-0.59(2) K. Assuming that the magnetic moments 
in the system are equal to 1/iB, we can estimate their 
concentration at the yttrium site (see discussion) as y = 
3 Ccw ffl 3 fce/16 /io Mbi where a is the lattice constant. This 
formula gives j/=0.60(l)%. 
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FIG. 2: (Color on-line) ZF (upper panel) and LF=50mT 
(lower panel) spectra as a function of temperature. The ini- 
tial muon beam polarization is parallel to the [100] crystal 
direction. The lines are fits as described in the text. 



the Z axis to be parallel to the initial muon spin di- 
rection, with the positron detectors centered along this 
axis. With this geometry the measured asymmetry (i.e. 
/iSR signal) is written as aoPz(t), where ao is the initial 
asymmetry (a characteristic of the spectrometer and the 
geometry of the sample), and Pz(t) is the longitudinal 
muon polarization function which contains information 
on the local magnetic fields in sample. In addition, three 
measurements were done with the Transverse Field (TF) 
geometry where B GXt is perpendicular to Z. 

It is important to know if there is muonium formation 
in Y 2 Ti 2 07, as this entity is known to diffuse coherently 
at low temperature in some systems. [24] Typical ZF and 
LF (50 mT) spectra measured at different temperatures 
are shown in Fig(2] The initial asymmetry is tempera- 
ture and field independent with a value ao ~ 0.25 typical 
for the GPS. This is a strong indication that there is no 
muonium formation in the whole temperature range. As 
expected for the absence of muonium, no re-polarization 
(a recovery of the missing asymmetry) by the LF is ob- 
served. The absence of muonium is further supported by 
the TF and ZF signals at 2.4 K (Fig. [3), where the initial 
asymmetry is found to be the same for the two spec- 
tra and the only frequency observed in TF is that of the 
applied field. This is further supported by Fourier Trans- 
form of the TF signal shown in the inset of Fig. [3] If an 
appreciable vacuum- muonium fraction was present, we 
would expect to observe an oscillation at a frequency « 
1037^5^/(271-) = 9.8 MHz ( 7m = 851.615 Mr/ads" 1 T" 1 ) 
in the Fourier spectrum. [24] and/or a reduction of the 
initial asymmetry due to the fast muonium precession. 
Therefore, we find no evidence of neutral muonium states 
in Y2Ti20y. The reason for the muonium absence in the 
pyrochlore oxides is not understood, and it is certainly a 
subject of much interest. 

We turn now to discussing the ZF and LF data in 



Fig. [2] where we observe a surprisingly high relaxation 
rate. A priori, electronic magnetism should not be 
present in Y2Ti 2 07 since both Y 3+ and Ti 4+ ions are 
non-magnetic. Therefore, only nuclear magnetic mo- 
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FIG. 3: (Color on-line) TF (0.745 mT) and ZF spectra mea- 
sured at 2.4 K. The inset shows the Fourier transform of the 
TF spectrum, where a peak is observed at 0.101 MHz corre- 
sponding to the applied TF. 

ments should be present in this system. However, as 
shown in Fig[l] the susceptibility shows a relatively 
strong paramagnetism. Also, as shown in Fig. [4j not 
even a field B eKt = 20 mT can completely decouple the 
relaxation of the /iSR signal (usually less than 3 mT 
are needed to decouple nuclear fields. [55] which we es- 
timated to be ss 20 /xT in Y 2 Ti207). Hence, we conclude 
that the relaxation is of electronic origin as proposed by 
Dunsiger.[2"T] 

To understand the origin of the electronic magnetic 
moments we first discuss the Y 2 Ti2C>7 crystal structure. 
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FIG. 4: (Color on-line) The field dependence of the LF data 
at 2.4K. The solid lies are fits as described in the text. 

In general, we note that the ternary pyrochlore com- 
pounds are of the general formula A 2 i? 2 07.[22] The B 
element can be a transition metal with a variable oxida- 
tion state. This gives the possibility of substitution on 
the B site (see Ref. and references therein). The 
structure also tolerates vacancies at the A and O sites to 
a certain extent. Around the A and B cations one finds 
8- and 6-fold coordination polyhedron of oxygen, respec- 
tively. Focusing on Y 2 Ti 2 07, in terms of the oxidation 
states we have Y 2 + Ti 2 + 07 + . Ti 4+ is diamagnetic since 
it is in a d° configuration. Depending on its coordination 
the Ti 4+ ionic radius is 74.4 or 88 pm, which is smaller 
than the Y 3+ ionic radius of 104 or 115.9 pm. Now we 
consider a magnetic configuration for at least one of the 
three involved elements. The obvious candidate is Ti 3+ 
which has about the same ionic radius as Ti 4+ . This 
ion is known to be rather unstable [27] and it has been 
suggested to locate it in the A site. [22] In this case we 
can write (Y^^Ti^T^C^, where we have neglected the 
substitution of Y 3+ in the B. This does not change our 
conclusions since Y 3+ is non-magnetic. [2"5l |2"9"] Ti 3+ has 
an electronic configuration d 1 and therefore, as a free ion, 
carries one Bohr magneton. Another possibility, which is 
less likely though cannot be rule out, is O vacancies in 
the structure. Such defects will also introduce two Ti 3+ 
magnetic moments for each vacancy. [28 Note that both 
types of defects will produce essentially the same effect 
on the implanted muons and cannot be distinguished in 
our measurements. 

We now consider qualitatively the ZF and 50 mT LF 
spectra of Fig. [2] In ZF and low temperatures we ob- 
serve a signal with a fast and a slow relaxing compo- 
nents. As the temperature is increased, the relaxation 
rate of the fast component decreases slightly and then 
increases until this component disappears above 150 K. 
The relaxation of the slow component decreases contin- 
uously with temperature. In the LF measurements the 
overall relaxation rate start increasing as the tempera- 



ture is increased, peaks at 100 K, and then it decreases 
to a small value above 150 K. Referring to the original 
work of Borghini et al. [30] and the latter works reviewed 
in Ref. [25], the temperature dependence of the spectra 
is typical for a muon trapping/detrapping effect. At low 
temperature the overall relaxation is small because a big 
fraction of muons are implanted "far" from Ti 3+ mag- 
netic defects. As the temperature is increased, the relax- 
ation increases because muons diffuse towards magnetic 
defects and become trapped. At even higher temper- 
atures, trapping/detrapping become faster leading to a 
decrease of the relaxation rate due to the fast fluctuating 
magnetic field sensed by the muons. The 50 mT field is 
sufficient to quench the relaxation except around 100 K, 
where a i? cxt = 200 mT was needed to suppress the relax- 
ation at this temperature. This is consistent with muons 
diffusing and accumulating near magnetic defects at this 
temperature, i.e. the trapping rate is higher than the 
detrapping rate. 

We base our analysis on a multi-state model and, for 
simplicity, we restrict ourselves to a two-state model. [30, 
31 In this model the muon can be diffusing in the undis- 
turbed regions or trapped near a magnetic impurity. The 
relaxation rate from muons in the trapped state is ex- 
pected to be higher since they are closer to magnetic 
defects and therefore experience stronger magnetic fields. 
Also, muons in undisturbed regions can diffuse and reach 
trapping sites (at a given trapping rate), while trapped 
muons can escape traps if the temperature is high enough 
(at a given detrapping rate). The muon polarization 
function for such a two-state model can be approximated 
by the sum of two relaxing signals. [31] 

P z (t) = /exp(-Ait) + (l-/)exp(-A 2 t), (1) 

where Ai and A 2 are the relaxation rates of the two 
components and / is the contribution of the first com- 
ponent to the full signal. At the limit of zero trap- 
ping/detrapping rates, Ai and A 2 represent the relax- 
ation rate of muons in traps and undisturbed regions, 
respectively, while / is the fraction of muons in traps. 
Eq. ([!]) provides a good qualitative description of the ZF 
data at all temperatures as shown by the solid lines in 
the top panel of Fig. [2] The temperature dependence 
of the three parameters in ZF is shown in Fig. [5] We 
find / ~ 0.20 and temperature independent below 40 K, 
indicating that the muon trapping and detrapping rates 
are constant below 40 K. Furthermore, the fact that for 
these temperatures the signal can be decoupled by small 
applied fields, is evidence that most of the muons can- 
not diffuse to reach a high field trapping site during the 
experimental time window (8.5/xs). Therefore, we con- 
clude that below 40 K the trapping and detrapping rates 
must be small. At higher temperatures / increases due 
to the enhanced muon trapping caused by a faster dif- 
fusion in the undisturbed regions, and then goes to zero 
at 135 K indicating that muons trap and detrap so fast 
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FIG. 5: The temperature dependence of Ai, A2 and / obtained 
from fits of the ZF spectra to Eq. (JTJ). f was set to zero for 
T>150 K. 



that they experience fast fluctuating magnetic fields (mo- 
tional narrowing) and relax following a single exponential 
behaviour. [3T] We want to point out that fitting the data 
with a temperature independent / does not produce sat- 
isfactory fits; and fitting with a power-exponential func- 
tion produces very bad results for temperatures above 
70K. Also, an implementation of the two state model as 
presented in Ref. [31] can not account for the observed 
temperature dependence of /; probably because of the 
assumption that the relaxation rate within each region 
(state) is exponential in the whole temperature range. 

To study the dynamic behavior at 2.4 K, we fit the LF 
data using the polarization function (see Fig. [4|, 



Pz(t) 



/ *P AJ (Ai,i/i, B cxt ,t) 

+ (1- f)P KT (A 2 ,v 2 ,B cxt ,t), 



(2) 



where P KT is the analytical approximation of a Gaussian 
Kubo-Toyabe function proposed by Keren. [32] By con- 
struction, this equation assumes two independent muon 
fractions and therefore it represents the two-state model 
in the limit of zero trapping and detrapping rates. In each 
region (or state) the muons sense random fields from a 
Gaussian distribution of width Aj. This field can fluc- 
tuate at a rate Vi due to fluctuations of the magnetic 
defects or hopping of the muon from one site to an other 
within the same region. In the fit, / was fixed to its 
ZF value (0.2) and B ext to the applied LF. The fit is 



good and completely captures the decoupling of the sig- 
nal, further supporting our assumption of small muon 
trapping and detrapping rates. The values of the fitted 
parameters are: A x = 3.00(5) mT, v\ — 2.4(2) /is -1 , 
A 2 = 0.318(5) mT and v 2 = 1.52(7) /xs" 1 . Ai, which 
correspond to the trapped state, is consistent with the 
dipolar field expected in a region of 8 A around a 
magnetic impurity (see discussion below) ; and v 2 imposes 
a maximum limit in the hopping rate of muons in the 
undisturbed regions. We want to note though that Vi 
and v 2 are very similar and, in fact, the data can be fit 
with a common fluctuations rate with no significant effect 
on the quality of the fit. Since the contribution from the 
paramagnetic defects to v\ and v 2 is the same, this indi- 
cates that the hopping rate in the undisturbed regions is 
probably much smaller than v 2 ■ 

One point that needs to be discussed is the relatively 
large / fraction found at low temperature. There are six- 
teen Y 3+ ions per cubic unit cell of volume V cc — a 3 . Let 
us denote y the percentage of magnetic defects relative 
to the Y 3+ population: y = x/2, i.e. there are 16y = 8x 
Ti 3+ ions in V^. Let us assume that a muon is trapped 
in a domain of relative weight / when implanted within a 
distance d from a defect, and also that the size and num- 
ber of these domains are small enough such that they 
do not overlap. Therefore, the volume around a defect 
in which a muon is trapped is = 8x(47r/3)d 3 . Since 
/ = V^/Vcc, we conclude that d = [(3/32ir)(f /x)]V 3 a. 
With / ~ 0.20 and y = 0.6% (see caption of Figjl|, we 
compute d = 0.79a = 8 A which is a reasonable value. [53"] 
This relatively long-range influence of a defect is due to 
the long range of both, the dipolar interaction and the 
nature of the elastic distortion field. [34] [35] 

In conclusion, we find that the muon hopping rate at 
2.4 K, in both stopping regions, is smaller than «2/xs _1 . 
This value is far below that observed for coherent diffu- 
sion in metals as well as the ~10 3 /is _1 proposed for coher- 
ent muon diffusion of positive-muons in Dv9Ti9O7.|20] 
Therefore we find it unlikely that coherent muon diffu- 
sion is present in Y2Ti207. Nevertheless, an appropriate 
theoretical microscopic calculation for the behaviour of 
positive-muons in this system is needed to completely 
rule out or confirm coherent muon diffusion. However, 
any theory should include the extended magnetic de- 
fects found by our measurements. In this regard, our 
observations add to the growing evidence that defects 
in pyrochlores are important to understand their mag- 
netic properties [13 l [28 ] I36H40] . and in particular those of 
Dy2Ti207.p9j Finally, an analysis to extract trapping 
and detrapping rates in the whole temperature range 
would have to follow the lines of that in Ref. [3TJ but 
allowing for non-exponential relaxation rates and/or go- 
ing beyond the two states model. 
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